Correlated dynamics of inclusions in a supported membrane 

Naomi Oppenheimer and Haim Dianian10 
Beverly and Raymond Sackler School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel 

(Dated: August 1, 2010) 

The hydrodynamic theory of heterogeneous fluid membranes is extended to the case of a mem- 
brane adjacent to a sohd substrate. We derive the couphng diffusion coefficients of pairs of membrane 
inclusions in the limit of large separation compared to the inclusion size. Two-dimensional com- 
pressive stresses in the membrane make the coupling coefficients decay asymptotically as 1 /r^ with 
interparticle distance r. For the common case, where the distance to the substrate is of sub-micron 
scale, we present expressions for the coupling between distant disklike inclusions, which are valid 
for arbitrary inclusion size. We calculate the effect of inclusions on the response of the membrane 
and the associated corrections to the coupling diffusion coefficients to leading order in the concen- 
tration of inclusions. While at short distances the response is modified as if the membrane were a 
two-dimensional suspension, the large-distance response is not renormalized by the inclusions. 
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I. INTRODUCTION 

Biological membranes are fluid bilayers made primarily of lipid molecules [ij . From a hydrodynamic point of view 
such a bilayer is a quasi-two-dimensional (quasi-2D) viscous liquid, whose molecules are constrained to flow along 
the 2D membrane surface while exchanging momentum not only among themselves but also with the surrounding 
three-dimensional (3D) solvent. Biomembranes contain also a high concentration of embedded inclusions — integral 
proteins and possibly also nanometer-scale domains — which perform key biological functions and are typically much 
larger than the lipids [H- Thus, from the same coarse-grained perspective, a biomembrane can be viewed as a quasi-2D 
suspension . We have recently used this perspective to investigate the correlated motion of proteins in a membrane 
freely floating in an unbounded liquid, and how the inclusions affect the response of such a membrane to stresses 
In many practical circumstances, however, the membrane is not free but attached to a solid substrate, such as the 
elastic scaffold of the cytoskeleton or an external surface to which a cell adheres. Additionally, in various experiments 
membranes are more easily studied when supported by a solid substrate. The aim of the current work, therefore, is 
to explore how the results of Ref. [3J are modified by the presence of such a nearby immobile surface. 

Two major approaches to the hydrodynamics of free membranes have been presented. The first, by SafFman and 
Delbriick (SD) [J| , models the membrane as a viscous liquid slab of width w and viscosity rjm/w, havingno-slip contacts 
at its bounding surfaces with two semi-infinite fluids of viscosity rji. The second, by Levine et al. [5|, considers the 
membrane as a vanishingly thin viscoelastic film embedded in an infinite viscous fluid. The dynamics of membranes 
embedded in a 3D fluid have been studied also using computer simulations The common key feature of these 

theories is the fact that the membrane does not conserve momentum in 2D, while the total momentum is conserved 
in 3D. Consequently, a length scale emerges, characterizing the crossover from a 2D-like membrane response, 
where stresses dominantly propagate through the membrane, to a 3D-like response, where the outer fluid governs the 
dynamics. This length is determined by the ratio between the 2D viscosity of the membrane and the 3D viscosity of 
the surrounding fluid For lipid bilayers is typically two to three orders of magnitude larger than the membrane 
thickness w, i.e., of micron scale. The work of Ref. Q is an extension of the SD theory to cases with more than one 
inclusion, where there are three lengths to consider: the lateral size (radius) a of the inclusion, the SD length k~^, and 
the distance r between the inclusions. The analysis was restricted to the limit a <C min(K~"'^, r), where complications 
related to specific details of the inclusions ^ are immaterial. We employ the same assumption in most of the current 
work as well. 

The introduction of an immobile surface breaks the translational symmetry in the directions parallel to the mem- 
brane and, hence, qualitatively changes the hydrodynamics of the system as its total momentum is no longer conserved 
Models appropriate for such a scenario can be divided into two groups. The first p^| - [l3 | adopts a phenomenolog- 
ical approach, describing the membrane as a 2D Brinkman fiuid jl5t] . i.e., introducing a term in the 2D hydrodynamic 
equation for the membrane, which leaks momentum at a certain fixed rate a^Vni (f^m ~ Vm/ Pm being the 2D kinematic 
viscosity of the membrane and p,n its 2D mass density). This sets a phenomenological screening length, a~^, beyond 
which momentum is lost to the substrate. In the second approach (16hil9i] the additional length scale is explicitly 
determined by the thickness /i of a fluid layer separating the membrane from the solid substrate. When fluid exists 
only between the membrane and the substrate (as in the case of a supported lipid monolayer) p^l - fisj . or when the 
membrane lies at the midplane between two substrates 0, H^, the large-distance effects converge to those of the 
phenomenological, Brinkman-like approach in the limit /i <C (i-e., for h much smaller than a micron). The 
momentum screening length thus obtained is the geometrical mean of the other two lengths, ~ (k~^/i)-'^/^ [17|. 

The current work extends the analysis of Refs. [16l - [l8| to the realistic scenario of a supported membrane with fluid 
on both sides, containing more than one inclusion. In Sec. |TT] we define the model and present the results for the 
hydrodynamics of a supported membrane, which will be useful for our analysis. Because of the three length scales 
in the model — k^^, h, and the interparticle distance r (assuming that a is much smaller than all three) — there 
are several asymptotic regimes to be considered, which are defined and discussed in Sec. |TT1 These results are used 
in Sec. IIIII to calculate the coupling diffusion coefficients of pairs of inclusions, which should be directly measurable 
in particle-tracking experiments. We proceed in Sees. IIVI and |V] to examine the effect of a finite concentration of 
inclusions on the response of the supported membrane to stresses (i.e., the effective viscosity of the membrane) and 
the resulting corrections to the coupling diffusion coefficients. In the practically useful limit of ft <^ we have 
been able to derive the large-separation coupling diffusion coefficients for inclusions of arbitrary size; these results are 
presented in Sec. IVIl The conclusions are summarized in Sec. I VIII 



II. MODEL 



Our model system is similar to that of Refs. |16l4l8| and is schematically depicted in Fig. [TJ A flat slab of viscous 
liquid (width w and viscosity rj-m/w) lies on the xy plane a distance h away from a flat rigid surface. The space 
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between the slab and the surface is filled with another fluid (viscosity r/f). Unlike the system treated in Refs. [l6l - [l8| . 
the space on the other side of the slab is occupied by a semi-infinite fluid of viscosity rjf as well. All fluids are assumed 
incompressible. The SD length is defined as 

K-' = ^. (1) 

In the slab are rigid inclusions of radius a. Although they are depicted as cylinders, their exact shape does not affect 
most of our analysis; in points where it does, this will be explicitly mentioned. We do not consider effects related to 
curvature and thermal fluctuations of the membrane la, |20I422| . 




FIG. 1. (Color online). Schematic view of the model system and its parametrization. 



Our goal is to characterize the response of the inclusion-decorated membrane to stresses and relate it to the coupled 
motions of two inclusions. A common way to characterize the membrane response as a fluid medium is through the 
velocity Green's function, G(r). This tensor gives the flow velocity v(r) of the membrane at the 2D position r due to 
a point force F exerted on the membrane at the origin in the xy plane, according to Vi{r) — Gij{r)Fj (with i, j — x,y 
and summation over the repeated index j). When the separation between two inclusions is much larger than their 
sizes (r ^ a), the pair mobility and pair diffusion coefficients associated with their coupled motions can be directly 
obtained from the velocity Green's function, as explained in detail in Sec. IIIII 

In unbounded 3D liquids G is the Oseen tensor We need the analogue of the Oseen tensor for the model 

system of Fig. [T] The velocity Green's function for a similar model of a supported monolayer, where the upper fluid 
is absent, was derived in Ref. [l6|. Its generalization to the case of two different fluids and two confining surfaces has 
been performed in Ref. [19]. We use this result while specializing to two identical outer fiuids, the upper one being 
semi-infinite. In Fourier space [G(q) = J (Pre^^'^''^G{r)] the Green's function is given by [1^ 

^"^^ " mq[coth{qh) + 2q/K + l] {^'' ~ ^) ' 
In the limit of an infinitely distant surface Eq. ^ coincides with the Green's function for a free membrane 

h->-qc _ 1 ( (5 - — '^^ 

In the other limit, of a vanishingly small h, it becomes 



G..'^G^,(q) = — Tf-T %-^). (3) 



Vm{q'+nq/2 + a^)[''' q^ 



with a = [n/ (2ft,)]^/^. The function in Eq. (|4]) has two poles at q± = n[-~l ± \/l — 8/(k/i)]/4, which in the hmit h ^ 
turn into q+ = ±ia. Hence, in this limit G coincides with the 2D Brinkman-like Green's function for an adsorbed 
membrane |10l - [l4| . 

7/m(g2 +a^) q- 



'^^ Gt, (q) = -7-f-5T - ^ ) ■ (5) 
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Thus, in the strongly adsorbed hmit the presence of the upper fluid merely adds a prefactor to the momentum 
screening length compared to the monolayer case studied earlier [13 [a-^ = (2«;-i/i)i/2 instead of {n-^hf/'^]. 

The Green's function of Eq. ([2|) serves as the starting point for the entire analysis to follow. It is beneficial, therefore, 
to begin by exploring the different asymptotic regimes that G(r) defines. Since it depends on three lengths — , 
h, and r — there are quite a few such regimes. Let us first rewrite it as a function of two length ratios. 



G 



9^J = 



{2n) 



1 



fc[cothfc+ (2K-V^)fc + 1] 



(6) 



thus making the asymptotics more transparent. 



A. Adsorbed regime: h <^ n ^ 

In this commonly encountered regime, where the distance to the substrate is much smaller than the SD length, the 
integral in Eq. ^ is governed for all values of r by small k, whereby coth k + {2k~^ /h)k ~ 1/fc + {2k~^ / h)k. Due to 
the argument presented below Eq. (|4]), the Green's function coincides with of Eq. ([5]). Inverting it to real space, 
one gets 



27r77n 



1 



27r r/mCt^r^ 



(7) 



where Kn are modified Bessel functions of the second kind. The first term in Eq. ([7]) is short-ranged, decaying 
exponentially with ar. It arises from shear stresses in the membrane, which get screened for r > a"^. The second 
term is long-ranged, decaying only as l/r^ and originating from compressive stresses in the membrane. The effect 
of such stresses on the steady fiow is that of an effective 2D mass dipole [2], whose magnitude is proportional to 
(f/mtt^)"^ h/r]f, independent of membrane viscosity. 

Thus, the adsorbed regime is subdivided into two regions reflecting different physics. In the adsorbed near region, 

response is governed by the yet-unscreened 2D shear stresses in the membrane. 
Expanding Eq. ([7]) in small ar, we have 



G^,ir) 



(r) 



1 



{- [ln(ar/2) + 7 + 1/2] <5„- + ^} + 0[(ar)Vr?„ 



(8) 



where 7 is Euler's constant. This result coincides with the Oseen tensor of a momentum-conserving 2D liquid, 
exhibiting the well known logarithmic behavior of this problem, with a cutoff length of a~^. In the adsorbed far 
region, r ^ = (2k~^/i)^/^ ^ h, the response is due to long-ranged compressive stresses, yielding 



G.,(r)^Gff(r) 



2m]fr'^ 



Si 



2rirj 



0[{ar) 



-1/2, 



(9) 



Note that, since in the adsorbed regime h < a ^ < k ^, the adsorbed near region includes distances r both smaller 
and larger than h, and the adsorbed far one includes distances both smaller and larger than k~^. 

B. Hovering regime: <C h 



In this regime, where the thickness h of the fluid layer between the membrane and the substrate is the larger length 
scale, the asymptotes of Eq. ^ depend on the value of r/h. For r <^ h the integral in Eq. © is governed by large 
k, whereby coth A: -I- {2k~^ /h)k ~ 1 -I- {2k~^ /h)k. The Green's function then coincides with that of a free membrane, 
of Eq. (|3]). This is the case studied in Ref. Q. Inverting G^ to real space d il yields 



G 



H,{nr) ~ _ 1 iY,{nr) - Y^inr)) + 

Kr 2 TT[Kry 



2Hi{Kr) 



nr 



Y2{nr) 



7r(fi;r)^ 



(10) 
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where are Bessel functions of the second kind, and i?„ Struve functions. 

The free behavior is subdivided into two regions having different physics. For r ^ <IC h we have the free near 
region, where Eq. (|10l) becomes 

G.,(r) ^ Gf;^(r) = [H>^r/2) + 7 + 1/2] % + ^} + 0{Kr/rj^). (11) 

As in the adsorbed near region [Eq. ([S])] the response in the free near region is governed by 2D shear stresses in the 
membrane. The difference is in the cutoff length, which in this case is rather than a"^. In the free far region, 
<^r ^h, Eq. (dU]) yields 

G.,(r) ^ Gf,(r) = + OK^r)-/^.]. (12) 

Both the 1 /r decay and the dependence on ?7f rather than indicate that the response in this region is due to shear 
stresses in the 3D fluid on the two sides of the membrane. 

In the last asymptotic region, <^ h r, although h is the larger characteristic length, the distance r is 
sufficiently large to make the fact that the membrane is supported rather than free come into play. In this supported 
region the integral in Eq. © is governed by small k, yet, unlike the adsorbed regime, cothfc + {2k^^ /h)k ~ l/fc + fc/3, 
independent of n. The resulting Green's function is 

which has two poles, both depending on h alone — i.e., h is the sole momentum screening length in this region. 
Inverting Eq. to real space and taking the limit r ^ h yields 

G,(r) ^ Gl^(r) = (s., ^ ^) + 0[hV{rj,r% (14) 

Thus, despite differences in the details (a~^ vs. h as the length scale for momentum loss), the asymptotic responses 
in the adsorbed far and supported regions [Eqs. ([9]) and (fT4|) . respectively] turn out to be identical, both arising from 
2D compressive stresses due to an effective mass dipole of magnitude ~ h/rji. Note, however, that the correction 
to the leading mass-dipole term in Eq. ([9]) is exponentially small, whereas in the supported region the correction is 
algebraic. This reflects the fact that in the adsorbed region the upper fluid is insignificant, whereas in the supported 
region it does play a role, albeit not a dominant one. 
The various spatial regimes are summarized in Table [H 



regime 


sub- region 


definition 


spatial dependence 


mechanism 


adsorbed 




h < k"^ 








near 


r<Q-i = (2k-i/i)1/2 


In(Qfr) 


2D shear 




far 


r > a^^ 


-2 

r 


2D compression 


hovering 




< h 








free near 


r < 


ln(«:r) 


2D shear 




free far 


K"^ < r < ft 




3D shear 




supported 


r > ft 




2D compression 



TABLE I. Summary of asymptotic spatial regimes and their corresponding notation. 



III. CORRELATED DIFFUSION 

The Green's function of Eq. (|6]) gives the membrane flow velocity at position r in response to a unit force exerted on 
it at the origin. In the limit r ^ a, addressed in this article, the same Green's function also gives the coupling mobility 
tensor, Bi2^ij{r) — i.e., the velocity of one particle due to a unit force acting on another, where the positions of the 
two particles are separated by the vector r, vi^i = i3i2,ij (r)-F2.j (with summation over the repeated index j). From the 
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mobility tensor the coupling diffusion tensor, Di2,ij(r), readily follows via the Einstein relation, Di2,ij = fcBT-Bi2,y , 
k^T being the thermal energy. The x axis can be defined, without loss of generality, along the line connecting the 
pair, r = rx. This choice leads, by symmetry, to Di2^xy = 0. The coupled diffusion of the two particles is then fully 
characterized by two coefficients: a longitudinal coupling diffusion coefficient, Di^{r) — Di2^xx{fx), and a transverse 
one, D^ir) — Di2^yy{rx). Thus, in summary, we have 



Di,{r > a) = kBTGxxirx), D^ir > a) = kBTGyy{rx) 



(15) 



The first coefficient is associated with the coupled Brownian motion of the pair along their connecting line, while the 
second — to the coupled motion perpendicular to that line. 



(16) 



where Ax^, Ay^ (/? = 1,2) are the displacements of particle /3 during time t. We shall now give the expressions for 
these coupling diffusion coefficients in the various asymptotic regimes (cf. Table 



A. Adsorbed regime: /i <C k ^ 



In the adsorbed regime we get from Eqs. ([7]) and the following coupling diffusion coefficients: 



^L(r) : 
DT{r) 



kBT 
27r?7m 
kBT 



Ki{ar) 



ar 



Koiar) + 



Ki{ar) 



This regime is subdivided into near and far regions. In the adsorbed near region Eqs. (m and ^ yield 



kBT 



whereas in the adsorbed far region we obtain from Eqs. ([9]) and ([TSl) 



[-In(ar/2) - 7± 1/2] , 



kBTh 
27r7]fr^ 



(17) 



(18) 



(19) 



In Eqs. (|18p and (IT51) the upper (lower) sign corresponds to the longitudinal (transverse) coefficient. Notice how at 
large distances, r ^ a~^, the coupling diffusion coefficients both decay as l/r^, suggesting the dominance of the 
aforementioned mass dipole. Momentum at such distances is lost to the solid substrate and the coupling is mediated 
by 2D compressive stresses in the membrane. 

The dependencies of the coupling diffusion coefficients on the separation between the inclusions for the adsorbed 
regime are shown in Fig. [5] along with their various asymptotes. 



B. Hovering regime: k <^ h 



The hovering regime is divided into the free (r <^ h) and supported (r ^ h) behaviors. In the free limit we get 
from Eqs. (HUl) and ^ 

< /i : Di^{r < /i) ~ \ Hiinr) - Yi(Kr) \ 

D^ir « /i) Uoi^r) - - 1 [Yo{nr) - Y2i>,r)] + ^] . (20) 

This behavior is further subdivided into near and far regions. In the free near region Eqs. (jlip and (jl5p yield the 
coupling diffusion coefficients as 

K'^^h: 7^LT(r<K~i)^-^[-ln(Kr/2)-7±l/2], (21) 
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FIG. 2. Longitudinal (A) and transverse (B) coupling diffusion coefficients as a function of interparticle distance for the 
adsorbed regime (h ^ k~^). Diffusion coefficients are scaled by ksT/rim and the distance by the momentum decay length 
. The full behavior [Eq. (I17p . solid] is shown together with those in the two asymptotic regions: adsorbed near region 
{Dl,t ~ Inr, dashed); adsorbed far region (Dl,t ±l/r■^ dotted). 



where, again, the upper (lower) sign corresponds to the longitudinal (transverse) coefficient. Note the similarity 
between this result and the one in the adsorbed near region [Eq. (|18p ]. In both cases the coupling is governed by the 
behavior of the membrane as a 2D fluid, the only difference being the cutoff of the logarithmic term. In the free far 
region Eqs. and ([T5|) lead to 

< /i : Di^{k^^ < r < /i) 



k^T _ kBTr]„ 



DTin « r « fe) - ^^^^^ ^ ^—212- (22) 



The coupling in this region is mediated by the outer 3D fluid, as reflected by the dependence of I?l on rjf and its 
spatial decay as 1/r. The transverse coefhcient decays faster (as 1/r^), since it arises from an effective force dipole 
proportional to ~ Vm/ilf [S]- This also leads to an unusual increasing dependence of Dt on membrane viscosity. 
All of the equations in the free limit, Eqs. (pHt - iP^ . coincide with those for a free membrane as derived in Ref. |3|. 
In the last region, the supported region, Eqs. (fT4|) and ([T5|) give the coupling diffusion coefficients 



Z?L.T(r » /i) ^ ±M1, (23) 

which are identical to those in the adsorbed far region [Eq. (flQl) ]. as they arise in both cases from the same physical 
mechanism (2D compressive stresses in the membrane). 

The dependencies of the coupling diffusion coefficients on the separation between the inclusions for the hovering 
regime, along with the various asymptotic regions, are shown in Fig. |31 



IV. EFFECTIVE RESPONSE 



The presence of inclusions in the membrane influences its response to stresses. In regular suspensions the response 
far from the point of perturbation is similar to that of the particle-free liquid but with a different prefactor, depending 
on the volume fraction of particles, (p. The modifled prefactor defines a modified, effective viscosity, ri'^^[(f)). 
For a 3D suspension of hard spheres the effective viscosity, to leading order in (/>, was calculated by Einstein [24] as 
rf^ = ri[l + (5/2)(/)]. The analogous calculation for a 2D suspension of hard disks yields [1^ rf^ = 77(1 + 20). As 
we have seen in Sec. |lll the inclusion-free supported membrane may exhibit a 2D-like response, a 3D-like one, or 
neither of the two. Hence, one expects a more complex modification of the response due to the presence of inclusions, 
depending on the various lengths in the problem. For sufficiently small distances momentum is conserved within the 
membrane, and we expect the response to be 2D-like. For sufficiently large distances transverse momentum is lost to 
the substrate, and the membrane is expected to behave neither as a 3D suspension nor as a 2D one. Heuristically, and 
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FIG. 3. Longitudinal (A) and transverse (B) coupling diffusion coeflicients as a function of interparticle distance for the hovering 
regime ^ h). Insets focus on the small-distance behavior. Diffusion coefficients are scaled by k-sT /rjm and the distance by 
the Saffman-Delbriick length . The full behavior [numerical integration of Eq. ((6| using nh = 10^, solid] is shown together 
with those in the various asymptotic regions: free near region (-Dl.t ~ Inr, dashed); free far region (Dl ~ l/^, Dt ~ 1/'''^, 
dash-dotted); supported region (Dl.t ~ ±l/r^ dotted). 



also based on the results for a free membrane f3|, we anticipate that the effective response of the supported membrane, 
having an area fraction (f) of disklike inclusions, G — > G°*^, will be obtained by the following transformation of the 
parameters: 



^Sf (0) 



r?m(l + 20) 
277f/77; 

a^a^^ic^) = [k"^/(2/i)] 



k(1 

1/2 



20) 
a(l 



(24) 



We now proceed to prove that this is indeed the case. 

Let us begin with an inclusion- free membrane and apply an in-plane, localized force density, F(5(r), at the origin. 
The resulting flow velocity field of the membrane is given by vf\Y) = Gij{r)Fj, where G(r) is the Green's function 
discussed in Sec. |IT1 Eq. ([5]). Next, let us consider the change in velocity at position r, (5v(r,r'), due to a single 
disklike inclusion located at r'. No force or torque are acting on the inclusion and, hence, its leading correction to the 
flow velocity is through the force dipole (stresslet) S that it introduces, Svi = Skj{v')dkGij{v — r'). There is a local 
relation between S(r') and the inclusion- free velocity field at r', given by 



Si 



27r?7n 



(0) 



(25) 



This membrane- analogue of Faxen's second relation [23] was derived in Ref. |3j] for a free membrane under the 
assumption na <C 1. It remains valid in the current case of a supported membrane, provided that a <C min(K~^, a~^). 

In the next stage we consider randomly distributed inclusions, occupying an area fraction (f> of the membrane. We 
restrict the calculation to the leading (linear) order in (f), where static correlations as well as hydrodynamic interactions 
between inclusions can be neglected. The average correction to the flow velocity is given then by integration over all 
possible positions r' of inclusions, multiplied by the uniform probability density of finding an inclusion centered at 
that position, (/)/(7ra^). 



j d\'Suj(v')duG,,{v-v') 



(26) 



The convolution in Eq. (|26p is conveniently handled in Fourier space. 



{5vM)) = ^5fe,(q)z(ZfcGy(q) = -20 



2g/K 



coth((7/i) + 2q/K+ 1 



Gy(q)F„ 



(27) 
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where in the last equation we have used Eqs. ([2]) and ([25)1 while neglecting the term of order (ga)^ in Eq. (|25l) . Writing 



~(o) 



(Svi) ~ G^fFj, we identify the effective Green's function as 



G 



off 



1 - 20 



2q/K 



coth{qh) + 2q/K + 1 



(28) 



It is readily verified that the same result is obtained from Eq. 1^ through the transformation defined in Eq. ([M)) and 
expansion to linear order in (j). 

The q-dependence of the renormalization factor in Eq. ((28|) reflects the aforementioned complex response of the 
inclusion-decorated membrane. In the g — > oo limit the prefactor becomes 1 — 2(j>, as in a 2D suspension '25'] . In the 
opposite limit of q — > the prefactor tends to unity — i.e., the inclusions have no effect on the large-distance velocity 
response of the membrane. To analyze the effective response in more detail, away from these two limits, one may 
invert Eq. (j28p back to real space in the desired limit. Alternatively, one can substitute the transformation defined 
in Eq. ((24|) in the various limiting expressions, already derived for G(r) in Sec. HI] [Eqs. (|7|)" ([T4|) ]. and expand those 
expressions to linear order in 0. Either of these two procedures yields the effective response of the membrane in the 
various regimes listed in Table HI We now examine the resulting expressions for the different regimes. 

In the adsorbed regime, h <C k^^, the substrate affects the response of the membrane. The correction to the Green's 
function of Eq. ([T]) due to the presence of inclusions is given in this regime by 



ft, < K 



^1. QCff 



G^ + SG'' 



6GUv) = 



27rry„ 



{~Ko(ar) + arKi(ar)) Sij — ariiTi (ar) 



(29) 



In the hovering free regions, k ^ <^ h and r <^ h, the membrane is insensitive to the presence of the substrate. The 
correction to the Green's function of Eq. ([T0|) for these free regions is 



< /i : G''*^ ~ G^ + 5G^ 
+ 



1 - {nrf 
4(Kr - 1) (nr)^ - 2 



7r(Kr)^ 



7r(Kr + 1) 

[H^iinr) + Yi{Kr)] - Ho{Kr) + Yo{Kr) 



r 



2 



(30) 



These two limiting cases (adsorbed and hovering free) are further subdivided into near and far regions. The near 
regions of both cases are governed by 2D shear stresses, leading to an effective response similar to that of a 2D 
particulate liquid with a leading logarithmic behavior. The differences arise from the different cutoff lengths — 
in the adsorbed near region and in the free near region. The different dependencies of these lengths on the 
membrane viscosity lead to slightly different concentration corrections as rj^ is modified to rf^ . In the adsorbed near 
region, r <C <C the correction to the Green's function of Eq. (|5]) is given by 

h < K"\r < a^^ : 



G 



cff 



(1 - 20)G: 



an 



— {[-(1 - 20) (ln(ar/2) + 7 + 1/2) + 0] % + (1 - 20)^} 



For the free near region, r ^ -C ft, the Green's function of Eq. pT|) is modified according to 



(31) 



r < K ^ < ft 



(1 - 20)G! 



27rr7,i 



{[-(1 - 20) {\n{Kr/2) + 7 + 1/2) + 20] S,, + (1 - 20)^} . 



(32) 



Equation p2|) has already been derived in Ref. Q for a free membrane. 

For all other regions we find no modification of the dominant term in the membrane response due to the presence 
of inclusions. This is because the response in these regions is governed by mechanisms which are unrelated to the 
propagation of 2D shear stresses in the membrane and, therefore, insensitive to membrane properties. In the adsorbed 
far region, ft <^ and r ^ a~^, and the supported region, <C ft <C r, shear stresses are lost to the substrate. 
The remaining compressive effect (effective mass dipole) is insensitive to the presence of inclusions, since the membrane 
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is assumed to be incompressible. Thus, the dominant response in these regions is given by the unperturbed Eqs. (|9]) 
and (O. In the free far region, <^r <^h, the response is dominated by 3D shear stresses in the adjacent fluid, 
which are obviously indifferent to the inclusions. Hence, in this region the dominant response remains equal to the 
unperturbed Eq. p^ . 

In the cases where the leading membrane response is not renormalized by the inclusions, there are nevertheless 
higher-order corrections which do depend on (j). To conclude this section we address these large-distance corrections. 
In the adsorbed far region the correction to Eq. © is exponentially small in ar and, therefore, negligible. In the free 
far region the leading correction to Eq. (IT2]) is obtained from Eq. ((30|) in the limit Kr 3> 1 as 



Thus, the correction is of higher order [1/r^ compared to l/r in Eq. (jl2p ] but still long-ranged. In the supported 
region we substitute in Eq. ((28)) the expansion coth(q/i) -I- 2q/ k ~ + 3qh and invert back to real space while 

assuming the limit r ^ h. This procedure yields the following correction to Eq. (|14|) : 



which decays algebraically but much faster than the unperturbed response (l/r^ vs. 1/ 



2\ 



V. CORRECTED PAIR-DIFFUSION COEFFICIENTS 



Substituting in Eq. ([T5|) the effective response functions, calculated in Sec. IIVI readily gives the corrections to the 
coupling diffusion coefficients to leading order in the area fraction of inclusions, 

Z^^^T - DuT + 5^L,T, (35) 

where Dl,t are the inclusion-free coefficients derived in Sec. IIIII We now provide the resulting expressions for (5Z?l,t 
in the various spatial regimes. 



A. Adsorbed regime: h <^ k ^ 

In the adsorbed regime we use SG'^ of Eq. in Eq. ([T5|) to get 

SDrir) ~ \Ko(ar) - arKi{ar)] . (36) 

These expressions are the corrections to the bare coefficients given in Eq. (1171) . In the adsorbed near region they 
reduce to 

h<^n-^: (5L»LT(r <a-i) ~<?;.-^[ln(ar/2) +7+1/2 T 1/2], (37) 

27r77,n 

where the upper (lower) sign corresponds to the longitudinal (transverse) coefficient. Equation ([37]) gives the correction 
to the bare coefficients of Eq. (fTS]). In the adsorbed far region, r ^ a~^, the corrections are exponentially small. 
More specifically, 5Dl ~ (ar)~^/^e~"'' and 6Dt ~ (ar)^/^e~"''. 

The behavior of the concentration corrections as a function of separation in the adsorbed regime is shown in Fig. H) 
Notice the range of positive correction to I?Tj i-e., the unusual increase in the transverse coupling due to the presence 
of inclusions. 



11 




ar 



FIG. 4. Corrections to the longitudinal (A) and transverse (B) coupling diffusion coefficients as a function of interparticle 
distance for the adsorbed regime [h <gi The corrections are scaled by (jfkBT/rjm and the distance by the momentum- 

screening length . The full behavior [Eq. H36|) . solid] is shown together with those in the two asymptotic regions: adsorbed 
near region ((5Dl,t ~ hir, dashed); adsorbed far region (5Dl,t ~ T(ar)=F'''^e-"", dotted). 



B. Hovering regime: k ^ <^h 

In the hovering free regions, ^ h and r <^ h, we substitute SG^ of Eq. ([50)1 in Eq. ([15]) to obtain 
K^^ < h : 



SDi^ir < h) : 
SDrir < h) 



-0 

2?7m 

,fcBr(i-M2) 

2?7n: 



Ho{kt) - + i {Y2{Kr) - YoiKr)) 

nr 2 



2^ r 



7rKr(l + Kr) 



(38) 



which are the corrections to the coefficients of Eq. ([20)1 . The hovering free behavior is further subdivided into near 
and far regions. In the free near region, the expressions in Eq. psp become 



which corrects Eq. ([^T|) . In the free far region Eq. (1551) 



27r?7in 
reduces to 



[ln(Kr/2) + 7 + 1 T 1/2] , 



,fcBT 1 



(39) 



(40) 



which are the corrections to Di^ and Dt of Eq. ([22]) . Note that in the free far region Dj- and SDt both decay as l/r^, 
whereas I?l has a slower decay (~ l/r) than its correction (~ 1/r^). Thus, the longitudinal coefficient in the free far 
region remains essentially unaffected by the inclusions. The results for the hovering free region coincide with those 
derived in Ref. 3] for a free membrane. 

In the last region, the supported region, where r ^ ft, 3> k^^, the dominant term in the large-distance response 
[Eq. ([Ti]) ] is insensitive to the properties of the membrane. The corrections to the coupling diffusion coefficients of 
Eq. ([23]), therefore, are of higher order. Substituting SG"" of Eq. ([34]) in Eq. ([15]), we get 



< h : 6Di^{r > ft) ~ -( 



12kBTh^ 

' 



48fcBrft3 



(41) 



The spatial behavior of the concentration corrections as a function of separation in the hovering regime is shown 
in Fig. [5] Notice again the broad range of positive correction to Dt^ where the transverse coupling increases with 
the concentration of inclusions. This is a consequence of the unusual dependence of the bare coefficient on membrane 
viscosity. [See Eq. ([^^ and the text below it.] 
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FIG. 5. Corrections to the longitudinal (A) and transverse (B) coupling diffusion coefficients as a function of interparticle 
distance for the hovering regime {k.^^ -C h). Insets focus on the small-distance behavior. Corrections are scaled by (jj^sT /rj^n 
and the distance by the Saffman-Delbriick length k"^. The full behavior [numerical inversion of Eq. (|28p using nh = 10^ solid] 
is shown together with the asymptotes for the free near region [SDi^^t ~ \nr, dashed) and free far region ((5Dl.t ~ Tl/*"^! 
dash-dotted). At much larger distances the supported region (JDl.t ~ dotted) sets in. 



VI. CORRELATED DIFFUSION OF LARGE INCLUSIONS 



The entire analysis so far has relied on the assumption that the inclusion size is much smaller than any other 
length in the system, a <C m.m.{K~^,h,r). Consequently, the coupling diffusion coefficients derived in Sec. IIIII were 
independent of the size and shape of the inclusions. In the adsorbed far region (r ^ 3> h) we can depart from this 
assumption and derive the large-separation coupling diffusion coefficients for large disklike inclusions. In principle, the 
calculation of pair mobilities for two large particles is technically hard — one needs to solve the flow equations with 
boundary conditions on the surfaces of the two particles as they move with prescribed velocities (or under prescribed 
forces). For the adsorbed far region, however, there is a scheme that bypasses this difficulty altogether. It is based 
on symmetry considerations and the knowledge of the exact flow field away from a single moving disk. 

In the adsorbed regime our model becomes equivalent to an effective 2D Brinkman fluid J^] — i.e., an incompressible 
fluid with momentum decay — as previously studied in Refs. [l3]-[i3- This is clearly seen in the corresponding velocity 
response function, Eq. ([5]), which contains a momentum decay length, a~^. 

First, let us recall, for the adsorbed regime, the far flow in the membrane due to a point force F. According to Eq. 
© this flow is the same as the one emanating from an effective 2D mass dipole of strength {h/r]i)F. Now consider 
the flow due to a force F applied to an isolated disk of radius ai, positioned at the origin. This problem was solved 
for a 2D Brinkman fluid in Ref. [ll|. Applying the result to our case, we get in the far field [r ^ max(ai,a~^)] the 
following dipolar flow: 



m{x) 



h 'm{aai) 



= 2 



xKo{x) + 2Ki{x) 
xKo{x)+4Ki{x)' 



(42) 



Thus, no matter how large the disk may be, the far flow remains equivalent to that induced by a 2D mass dipole; the 
only dependence on the particle size is through the dimensionless prefactor m{aa). 

In the next step we identify the tensor multiplying F in Eq. with the coupling mobility of two different particles 
— one of radius ai at the origin and another of vanishing radius, 02 — t- 0, at r. Due to the symmetry of the coupling 
mobility, the same tensor also gives the velocity of a disk of radius ai, positioned at r, due to a point force F applied 
to the membrane at the origin. Yet, the latter is the velocity acquired by a particle of radius oi as it is embedded in 
a flow caused by a mass dipole of strength {h/r]f)F. All we need to do now to get the velocity of a disk of radius ai 
due to a force F exerted on a sufficiently distant disk of radius 02 is to increase the mass-dipole strength at the origin 
from {h/ri{)F up to {h / rji)m{aa2)'P , the effective mass dipole created by F when it is applied to a particle of radius 



13 



a2- Hence, the large-distance coupling mobility is given by 



h m{aai)m{aa2) 2r,r^ , 
-«i2/y(rj = 7^-3 \°i3 3- ■ (.4^) 



As explained in Sec. IIIIl the coupling diffusion coefficients can be readily obtained from Eq. (|43| as Di^{r) — 
kBTBi2^xx{rx) and £>t(?') = kBTBi2^yy{rx), yielding 

1 N „ . N k'QT'm{aai)m(aa2)h , , 

r»max a-\ai,a2 : i^L.xW ^ ±^ ^ , ' 44 

where the plus (minus) sign corresponds to the longitudinal (transverse) coefficient, and m{x) is defined in Eq. (j42p . 
Equation (j44p gives the large-distance coupling diffusion coefficients in the adsorbed regime for two disklike inclusions 
of arbitrary radii. 

In the limit of small inclusions, aap <^ 1 {(3 = 1, 2), we have m(Q!a^) ~ 1, and the result of Sec. IIIIl [Eq. is 
recovered. In the opposite limit of large inclusions, m{aap) ~ 2, and, therefore, 

r>afl>a-i: Di^T{r) - ± " „ . (45) 

7r?7fr^ 

Thus, going from small to large inclusions changes the large-distance coupling coefficients by a mere factor of 4. 
Interestingly, the results for very large inclusions are again independent of particle size and shape. We return to this 
surprising finding in the next section. 



VII. CONCLUSIONS 



The aim of this work has been to characterize supported membranes as effective heterogeneous fluids. The existence 
of several length scales in the problem leads to various regimes that are governed by different physical mechanisms 
and exhibit different effective dimensionality (either 2D or 3D); see Table HI We have provided predictions for the 
coupling diffusion coefficients of inclusion pairs in those various regimes, as well as their leading dependence on the 
concentration of inclusions in the membrane. These predictions can be directly checked in two-point microrheology 
experiments using Eq. (jl6p . 

Since the SD length, k^^, is typically of micron scale, common supported membranes should belong in the adsorbed 
regime, h <C k~^, which is treated in Sees. IIII XI and IVXl Moreover, the limit of small inclusion size, assumed in those 
sections, should be generally valid, since the requirement is that a be much smaller than ~ (k~^/i)^/^ rather 
than the stricter condition a <^ h. Hence, we expect this limit to hold for common membrane inclusions even in cases 
where the distance to the substrate is of the order of the inclusion size (say, a few nm only). In this common scenario 
of a <C ^ the substrate is predicted to strongly suppress the large-distance correlations as compared to a 
free membrane. Comparing the results of Sec. IIII Al with those of Ref. 3] (or with the equivalent results for the free 
far region in Sec. IIII B|) . we find suppression of the longitudinal and transverse coefficients by factors of order h/r and 
Kh, respectively. Nonetheless, the correlations always remain long-ranged, their fastest possible decay being as 1/r^. 

For nm-scale separation between membrane and substrate, which is comparable to the membrane thickness, the 
substrate may affect the membrane properties. As long as the membrane remains fluid, such interactions are expected 
to merely modify the effective membrane viscosity, and the theory presented here should remain valid. A more serious 
concern is the possible breakdown of the bilayer description as a uniform slab, which is inherent in the Saffman- 
Delbriick model and the current work. At sufficient proximity to the surface the dynamics of the two membrane 
leaflets might decouple. This will occur when the friction between the lower leaflet and the solid surface exceeds the 
one between the two leaflets. The characteristic coefficients for these two competing drags are, respectively, rjf/h ^ 10^ 
N-s/m'^ (for h ^ 1 nm) and 10® N-s/m'^ ^26*1. Thus, for all relevant separations h, relative motion of the leaflets should 
not play a significant role, and the bilayer can be considered as a single fluid medium. 

There may be cases where the inclusion size is comparable to or larger than a^^ — for example, when the inclusion 
is a colloid particle or a membrane domain. We have presented expressions for the large-distance coupling coefficients 
in this case as well (Sec. IVI[) . These results have been derived for the specific case of disklike inclusions, yet in both 
limits of small and large aa they become independent of the size and shape of the inclusions. The origin of this 
surprising universality is that, in the adsorbed regime, the membrane responds to any size and shape of perturbation 
sufficiently far away, as if the perturbation were a mass dipole. The 'effective inclusion' — i.e., the region around the 
perturbation whose dynamics determines the strength of that mass dipole — is limited in both cases of very small 
and very large a by the momentum decay length a~^. 
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Concerning the effective response of tlie supported membrane as a function of the area fraction of inclusions, we 
have found that the membrane viscosity is modified according to the law for 2D suspensions [25j, rf^ = 77m (1 + 20), 
yet this modification should be included also in the parameters k and a [Eq. (|24|) ]. The combined effect is that 
there is no renormalization of the large-distance membrane response with increasing (j). The underlying physics is 
that transverse momentum is not transferred through the membrane over large distances and, hence, the response is 
insensitive to changes in the membrane viscosity. This insensitivity holds already for distances much smaller than the 
largest length scale in the problem at hand. For the common, adsorbed regime, h <^ k~^, it is valid for r 3> a^^, 
as momentum is first lost to the substrate. For the hovering regime, ^ ft., the dependence on (f> disappears for 
r 3> as the propagation of stresses becomes dominated by the outer fiuid. 

The subtle effect of increasing the area fraction of inclusions on the membrane viscosity and the governing length 
scales influences also the corrections to the coupling diffusion coefficients (Sec.|V|. At sufficiently short distances the 
leading logarithmic terms are corrected as if the membrane were a 2D suspension, yet at large distances all corrections 
to the coupling coefficients vanish. 

The validity of the theory presented here is limited in several important respects. Our results concerning the 
coupling diffusion coefficients are all valid only for large separations, r ^ a. In addition, the coupling coefficients 
derived in Sec. I VII for large inclusions apply only in the adsorbed far region, r 3> a^^ 3> h; the corresponding 
expressions for large inclusions in the other asymptotic regions are unknown. We have restricted the analysis of 
the effective response and corrections to the coupling coefficients to the leading linear order in the area fraction of 
inclusions, <j>. Deviations from the theory are expected, therefore, as cj) becomes appreciable. Nevertheless, as is 
clear from the discussion above, our main qualitative results — in particular, the insensitivity of the large-distance 
response to (f) — are expected to be valid for all values of 0, so long as the membrane remains fluid. Finally, we have 
not considered membrane fluctuations, which may have a subtle interplay with the diffusion of membrane inclusions 
[1, [131 . As membrane fluctuations are suppressed by the presence of a nearby surface [13, [2l| , we do not expect our 
main results to be significantly influenced by such effects. 
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